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This talk is an attempt to summarize some of the first results obtained at RHIC. I 
discuss the significance of these measurements for establishing the properties of hot and 
dense QCD matter and for understanding the dynamics of the theory at the high parton 
density, strong color field frontier. 

1. Introduction 

1.1. What is RHIC? 

Relativistic Heavy Ion Collider (RHIC) at Brookhaven National Laboratory on Long 
Island, New York, began operation in 2000 culminating over ten years of development and 
construction, and a much longer period of theoretical speculations about the properties 
of hot QCD matter produced in nuclear collisions in the collider regime. 

RHIC’s 2.4-mile rings contain superconducting magnets, which operate at minus 451.6 
degrees Fahrenheit, 4.5 degrees above the absolute zero. RHIC collides two intersecting 
heavy ion beams at center-of-mass energy of up to 200 GeV/A (at luminosity of up to 
10 26 sec _1 cm 2 , which can be further increased in the future), and polarized proton beams 
at c.m.s. energy of up to 500 GeV. The total energy in the gold-gold collision thus reaches 
40 TeV, which is at present the World’s record collision energy. In the pp mode, the unique 
possibility offered by RHIC for the first time is the study of double spin asymmetries and 
other spin observables. 

1.2. RHIC experiments 

There are currently five experiments at RHIC: BRAHMS, PHENIX, PHOBOS, STAR, 
which can operate in both heavy ion and pp mode, and PP2PP, which aims exclusively 
at the study of elastic and diffractive pp interactions. These experiments are different in 
their focus; however all of them also provide the information about the basic properties 
of the collisions. The data from different collaborations can thus be cross-correlated and 
independently verified. 

*Invited talk given at the XX International Symposium on Lattice Field Theory “Lattice 2002”, 
Cambridge, MA, June 24 - 29, 2002. 
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The BRAHMS experiment is designed to measure charged hadrons over a wide range 
of rapidity and transverse momentum. PHENIX is a large versatile experiment designed 
to study the production of charged and neutral hadrons, leptons, and photons. The 
PHOBOS experiment is aimed at measuring the global properties of nuclear collisions, 
and at the study of fluctuations and correlations in the production of hadrons. STAR, 
one of the two large-scale RHIC experiments, focuses on measurements of the production 
of various hadron species over a large solid angle and is well suited, in particular, for the 
study of multi-particle correlations and fluctuations on the event-by-event basis. 

1.3. RHIC physics: QCD 

RHIC is a machine entirely dedicated to the study of Quantum Chromo-Dynamics - 
the theory of strong interactions. This includes the study of matter at unprecedented 
energy densities and understanding the structure of the proton, in particular, its spin 
composition. Since the spin program at RHIC at the time of this talk is at the very 
beginning, I will concentrate in this talk entirely on the heavy ion program. 

Asymptotic freedom of QCD IB- IB ensures that the dynamics of quarks and gluons 
at sufficiently high density can be addressed by weak coupling methods. This includes 
both the thermalized quark-gluon plasma at high temperature and the wave functions 
of the colliding nuclei described, at small Bjorken x, by parton saturation and the Color 
Glass Condensate [ |3], [4], |5], ||. 

2. Looking at the first RHIC data 

2.1. Global observables 

Global observables are the most general characteristics of the collision, including par¬ 
ticle multiplicity, its dependence on the centrality of the collision and on rapidity, and 
azimuthal distribution of the produced particles. The centrality is determined by the im¬ 
pact parameter in the collision; since this is not a quantity which can be measured directly, 
centrality is usually determined in terms of a certain cut in the multiplicity distribution; 
e.g., a 0 — 10% centrality cut means that out of a given sample, 10% of the events which 
have the highest multiplicity have been selected. It is convenient to characterize central¬ 
ity in terms of the number of “participants” - the nucleons which underwent an inelastic 
interaction in a given collision. Glauber theory [ 0 can be used to correlate a certain 
centrality cut with an average number of participants (for an explicit set of formulae for 
nuclear collisions, see e.g. [ ||, [0). This procedure can be independently verified by 
measuring the energy carried forward by spectator neutrons; for that purpose all RHIC 
experiments are equipped by Zero Degree Calorimeters. 

2.1.1. Multiplicity 

Multiplicity in heavy ion collision tells us which fraction of the collision energy is in- 
elastically transferred to secondary particles. Theoretical expectations for hadron multi¬ 
plicities at RHIC varied by factor of five, and the experimental verdict was thus eagerly 
awaited. After the commissioning of the machine, the first multiplicity results did not 
take long to come; they are shown on Fig. [l] in comparison to the multiplicity previously 
measured in pp and pp collisions. The measured multiplicity appeared much smaller than 
most theoretical predictions. Is this disappointing? To answer this question, let us recall 
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Figure 1. Charged multiplicity per participant pair in central Au Au collisions as measured 
at RHIC is compared to the multiplicity in pp and pp collisions; from 
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that, by the very definition, an incoherent superposition of independent nucleon-nucleon 
collisions yields multiplicity equal to the number of collisions times the multiplicity in 
NN collision. This trivial statement holds also in the presence of elastic rescatterings. 
Indeed, according to so-called AGK cutting rules 
nuclear cross section is given by 


TOU of multiple scattering theory, the 
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where the nuclear overlap function is 


( 1 ) 
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and the nuclear thickness function T A (b ) = dzp A {b , z) is the integral over the nuclear 

density. Integration over impact parameter b yields 
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and correspondingly the particle multiplicity would scale as 
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Using the numbers of collisions (~ 1050) and participants (~ 340) in central (0 — 6% 
centrality cut) Au Au collisions from Glauber model calculations [ ;9|, we would thus 
conclude that AuAu multiplicity per participant pair should exceed NN multiplicity by 
















4 


factor of 3. Instead, the data at the highest RHIC energy of y/s = 200 GeV show the 
difference of only about 50%. Given that any inelastic rescatterings in the final state 
can only increase the multiplicity, we therefore have an experimental proof of a high 
degree of coherence in multi-particle production in nuclear collisions at RHIC energies. 
The diagrams which allow to evade the AGK theorem are Gribov’s “inelastic shadowing” 
corrections [ P"]J which correspond to the excitation of high-mass intermediate states in 
multiple scattering; the process thus no longer can be decomposed in terms of elementary 
NN interactions. In parton language, these contributions correspond to multi-parton 
coherent interactions. 
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Figure 2. Centrality dependence of the charged particle multiplicity near mid-rapidity in 
Au + Au collisions at A fs = 130 GeV and 200 GeV, from [ [HJ. 


2.1.2. Centrality dependence 

The dependence of multiplicity upon the number of participants discussed above can 
be established by selecting different centrality cuts. The result is shown in Fig.|2|; one can 
see that the multiplicity per participant pair increases with centrality, but not quite as 
fast as it would if the NN collisions were independent. 

If we decompose the multiplicity measured in NN collisions at some energy a fs into 
a fraction X(s) coming from “hard” processes, and the remaining fraction 1 — X(s) 
coming from “soft” processes, and assume that in nuclear collisions “hard” processes are 
incoherent and thus scale with the number of collisions, whereas “soft” processes scale 











5 


with the number of participants 

[0 


14fl, we arrive at the following simple parameterization 


dnAA = [(1 - X(s)) (1 V*) + X(s) (JV„„)] dnNN 


drj 


di] 


( 5 ) 


which describes the data quite well. In the framework of perturbative QCD approach, one 
has to assume that the coefficient X (s) is proportional to the mini-jet production cross 
section, and thus grows with energy reflecting the growth of the parton distributions at 
small Bjorken x, X(s ) ~ [: xG(x )] 2 , with x ~ 1/a/s. Therefore one expects [ [L5|] that the 
centrality dependence should become increasingly steep as the yfs increases (for the latest 
development, see however [|l(J). This increase is not seen in Fig.|], which in the lower panel 
shows that the ratio of the distributions at a/s = 200 GeV and a/s = 130 GeV is constant 
within error bars. The almost constant ratio appears consistent with the prediction [ [17 


Q based on the ideas of parton saturation, where the increase of multiplicity stems from 
the running of the QCD coupling constant determining the occupation number ~ 1 /a s 
of gluons in the classical held. The forthcoming results at a/s = 20 GeV will further 
constrain the production dynamics. 


2.1.3. Rapidity distributions 

Distributions of the produced particles in the emission angle 6 (with respect to the colli¬ 
sion axis), or pseudo-rapidity 77 = — ln[tan( 6 l / 2 )] provide another important characteristic 
of the collision process. Two features of RHIC results (see Fig. |3j) are especially notewor¬ 
thy: i) the distributions do not exhibit scaling in 77 ; ii) the deviation from NN results is 
maximal in the central rapidity region whereas the shapes of the AA and NN distributions 
are similar in the fragmentation region. Moreover, when corrected for the different beam 
rapidity 77 —► rj — ybeam , distributions at different energies exhibit approximate scaling in 
the fragmentation region (“limiting fragmentation”). 

2.1.4. Azimuthal distributions 

Of great interest and importance are the distributions of the produced hadrons in the 
azimuthal angle. Indeed, if all of the NN collisions were independent, there would be no 
reason to expect asymmetry in the distribution of the produced hadrons in the azimuthal 
angle (measured with respect to the reaction plane). This is why the azimuthal asymmetry 
represents a sensitive test of the collective effects in nuclear collisions. The azimuthal angle 
distributions of the produced hadrons are usually expanded in harmonics in the following 
wayQ: 

dN 

-= 1 + 2v\ cos p + 2v 2 cos 2ip + ... (6) 

dip 

Fig. [| shows the extracted from RHIC data coefficient v 2 {y 2 / 0 in the language of the 
held corresponds to “elliptic how”). One can see that the asymmetry of the azimuthal 
distribution is quite sizable, and for peripheral collisions (small multiplicity n c h/n max ) 
reaches about 35%. 

3 The absence of the terms proportional to sinny is the consequence of parity conservation; it would be 
interesting to search for their presence in the data in view of the speculative scenarios allowing for P and 
CP violation in hot QCD [ |2o| . 
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Figure 3. Pseudo-rapidity distributions of charged particles from Au + Au collisions at 

18| . Solid line is the prediction based on parton 
T7tl , and dashed line is the multi-phase transport model calculation [ |Hj . The 


y/s = 200 GeV (open circles), from 
saturation 


pp distribution rescaled by (N part / 2) is shown by stars. 


This effect certainly indicates the presence of collectivity in nuclear collisions, and comes 
out quite naturally in hydro dynamical calculations which assume complete thermalization 
in the final state [ |22[, [ |23|. However, final state effects are not the only possible origin 


of the azimuthal asymmetry; indeed, as we have discussed above, the high degree of 
coherence in the initial state signaled by the multiplicity measurements, in the parton 
saturation scenario, introduces strong correlation between the transverse momentum of 
the parton and its transverse coordinate in the wave functions of the nuclei. When 
the nuclei collide, this effect mimics at least a part of the asymmetry usually ascribed 


exclusively to final-state interactions 


25| . The magnitude of the elliptic flow which 


can originate from the coherence in the initial state is still a subject of ongoing research 
and debates. 

The dependence of the elliptic flow on the transverse momentum of the hadrons presents 
a puzzle [ |2T[ ; the value of v 2 is seen to first increase with p t , and then saturate. This 
contradicts to hydrodynamics, predicting a monotonic increase of V 2 with p t ; of course 
hydrodynamics cannot be trusted at high p t anyway because the density of hard particles 
is too small to allow for a meaningful statistical description. Energy loss of the produced 
jets can contribute to the azimuthal asymmetry at high p t (see, e.g., [0]), even though it 
is not yet clear if the magnitude of the effect can be reproduced under realistic assumptions 
about the density of the medium and the jet interaction cross section [ ET . 


2.1.5. Hadron abundances 

The measurements of yields of different hadrons at RHIC hold many surprises. Of 
particular interest is the fact that even at RHIC energies the asymmetry between baryons 
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Figure 4. Azimuthal anisotropy of hadron production in Au + Au collisions at at y/s = 130 
GeV; i >2 is the weight of the second harmonic, cos 2 p, in the particle distribution in the 
azimuthal angle ip ; from 
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and antibaryons is still sizable, with p/p ratio about 0.65 
of baryon number to quite small values of Bjorken x ~ ICE 2 . 


281. This signals the diffusion 


2.2. Hard processes 


2.2.1. Suppression of high p t particles 

Jet energy loss was among the first signatures suggested for the diagnostics of the hot 
quark-gluon matter [ |29], [3(| [H]]. This is why the measurements of the high p t hadron 
production excited a lot of interest. Indeed, the experimental results are striking - as can 
be seen in Fig.[5], the yield of high p t hadrons is drastically reduced with respect to what 
is expected for incoherent production in NN collisions. This behavior is very different 
from what was previously observed in Pb — Pb collisions at CERN SPS and in a — a 
collisions at CERN ISR (see Fig.^j). Does this important discovery signal jet energy loss 
in the quark-gluon plasma? The answer to this question can be given after we know 
the results of the forthcoming measurements in p(d)A collisions, which would allow to 
separate clearly the effects coming from the initial state. 


2.2.2. B/ 7T puzzle 

Another striking puzzle at RHIC is a rapid increase of the baryon-to-pion ratio in 
central Au — Au collisions at high p t [ The growth of this ratio is expected in 

the hydrodynamical scenario, in which equal velocity of the expanding parton “fluid” 
implies higher transverse momentum for more massive particles. However, the validity 
of hydrodynamical description is dubious at high p t where the density of particles is too 
small. If we assume that minijet fragmentation is the leading production mechanism of 
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Figure 5. The ratio of transverse momentum distributions of charged hadrons and neutral 
pions in Au + Au and pp collisions at -y/s = 130 GeV; from [0- 


high p t particles, then the growth of the B/ti ratio implies that minijet fragmentation is 
severely affected by the medium. Another scenario [ [34J attempts to explain both B/i r 
puzzle and a large value of baryon asymmetry B/B ^ 1, and invokes the contribution of 
non-perturbative gluonic junctions in nuclear collisions [ 

2.2.3. Charm production 

The production of heavy flavors and quarkonia represents a very important and exciting 
part of RHIC program. While these studies will benefit in the future from increased 
luminosity and improvements in the detectors (allowing, in particular, to reconstruct the 
decay vertex of heavy hadrons), the first measurement of charm production cross section 
has been already reported [ [I6[]. This has been done by deciphering the charm decay 
contribution to the single electron spectrum - see Fig|| Of particular interest is the 
fact that while the production cross sections of light hadrons, as discussed above, show 
strong nuclear effects, the cross section of charm production, within the error bars of the 
measurement, scales with the number of AAV collisions. These results may imply much 
smaller, in comparison to light quarks, energy loss of heavy quarks [ |37] stemming from 
the suppression of the gluon radiation at small angles (“dead cone effect”). 
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Figure 6. The background-subtracted electron spectra for minimum bias (0 — 92%) and 
central (0 — 10%) collisions compared with the expected contributions from open charm 
decays; from [ |36| . 


3. What have we learned so far? 

It is too early to assess the implications of RHIC results; however, it is becoming 
increasingly clear that they challenge most, if not all, of the existing theoretical dogmas. 
A coherent and convincing theory describing all of the observed phenomena directly in 
terms of QCD still has to be born. However, we can already conclude that many of the 
observed phenomena clearly manifest collective behavior; nuclear collisions at RHIC are 
not an incoherent superposition of nucleon-nucleon collisions. 

The measured particle multiplicities and transverse momentum spectra allow to esti¬ 
mate initial energy density at the early moments of the collision; a typical value inferred 
in this way is about 20 GeV/fm 3 4 5 6 (see, e.g., [ H). The dynamics of strongly interacting 
matter at such energy density (exceeding the energy density in a nucleus by over two or¬ 
ders of magnitude!) should be described in terms of quarks and gluons, and the collective 
phenomena observed at RHIC thus directly reflect the properties of high density QCD. 
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